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Abstract 

We study the transport properties of two electrons in a quasi one- 
dimensional disordered wire. The electrons are subject to both, a 
disorder potential and a short range two-body interaction. Using 
the approach developed by Iida et al. [ Ann. Phys. (N.Y.) 200 
(1990) 219 ], the supersymmetry technique, and a suitable truncation 
of Hilbert space, we work out the two-point correlation function in the 
framework of a non-linear a model. We study the loop corrections 
to arbitrary order. We obtain a remarkably simple and physically 
transparent expression for the change of the localization length caused 
by the two-body interaction. 

PACS numbers: 72.15Rn 71.30+h ] 

Non-interacting electrons in a one-dimensional random potential are lo- 
calized. The presence of a short-range two-body interaction (a screened 
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Coulomb potential) affects the localization length. For two interacting elec- 
trons in one dimension, this fact has been clearly established in a series of 
numerical and analytical investigations (see Refs. and references 

therein). The localization length increases independently of the sign of the 
interaction when the electrons move together at short distance 0. 

In the present work we address the more realistic and more demanding 
case of two electrons in a quasi one-dimensional disordered wire interact- 
ing via a short-range interaction. Using an analytic method developed re- 
cently f§], a reduction of Hilbert space and supersymmetry, we achieve a 
complete analytical solution of the problem. We show that the two-body 
interaction affects the localization length via the level density. We establish 
a criterion for the onset of interaction-induced effects on the localization 
length. Our final analytical expression allows for the numerical determina- 
tion of the localization length for any given two-body interaction. To the 
best of our knowledge, this is the first time that an analytical solution of a 
problem involving both, disorder and interaction, has been achieved. 

Model. We consider a quasi one-dimensional wire of length L which 
consists H of K slices labelled a,b,c, . . . , K. The surfaces connecting neigh- 
boring slices are transverse to the axis of the wire. Eventually we take the 
limit where K — > oo and where the longitudinal extension of the slices goes to 
zero. In each slice we use an arbitrary basis of single-particle states labelled 
\aj) with j = 1,2, ... ,1 and Fermionic creation and annihilation operators 
a' a j and a a j, respectively. We later take the limit I — > oo. 

The Hamiltonian H is the sum of three terms, H = H + Hi + H 2 . In 
every slice, disorder is simulated by an ensemble of random single-particle 
Hamiltonians belonging to the unitary ensemble, 

Ho = J2 h i? a li a w ■ (!) 

aij 

The complex random variables h°j have a Gaussian distribution with mean 
value zero and the following non-vanishing second moments, 



hifk^fj — —Saa'Sii'Sjj' . (2) 

The overbar denotes the ensemble average and A has the dimension of en- 
ergy. The disorder Hamiltonians in different slices are uncorrelated. Hopping 
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between neighboring slices is described by the hopping term 

H i = v l a li a a+u + h.c.) . (3) 

ai 

The short-range two-body interaction acts only between electrons in the 
same slice and reads 

H 2 = Y1 H WijiiyaliOtljaiajiaaii . (4) 

a i<j,i'<j' 

The two-body matrix elements are antisymmetric in the pairs (ij) and (i'j') 
and Hermitean. Both v and the u^y's are non-random quantities indepen- 
dent of the slice label a. For the single-electron problem, it has been shown 
via supersymmetry |J that H + Hi defines the same non-linear sigma model 
as obtained from a continuum model using a kinetic term and a random po- 
tential. This is why our model is generic. 

The Hilbert space of the two-electron problem is spanned by the or- 
thonormal states \abfi) where a < b and where /i stands for the pair (i, j) 
with i < j for a = b. For a < b fixed, the number N a b of states is I 2 while 
for a = b, N aa = 1(1 — l)/2. The totality of states \abfi) with \i — 1, . . . , N a b 
is referred to as the box (ab). The effect of the two-body interaction on 
localization properties in the one-dimensional case is strongest when the two 
electrons stay close to each other [j3|. This fact, and the need to reduce the 
dimension of Hilbert space, cause us to consider only the IK — 1 boxes where 
the two electrons are either located in the same slice a (this is box (aa)), or 
in two adjacent slices a and a + 1 (this is box (a(a + 1))), with a = 1, . . . , K. 

Supersymmetry. Localization properties can be read off the two-point 
correlation function 

C(n) = \ (a(a + l)fi\(E+ - H)~ l \ (a + n)(a + n + l)v) | 2 . (5) 

To be independent of edge effects, we choose 1 -C a (a + n) K. 
For large n, C(n) should decay exponentially in n. In proper units, the 
coefficient in the exponent defines the localization length. We calculate C (n) 
using supersymmetry || [/]]. 

Because of the Gaussian distribution of the coefficients h^f , the ensemble 
average in Eq. (|) is entirely determined by the second moments of the matrix 
elements (abfi\H \aba) of H with b = a, (a + 1), i.e., by the quantities 

•^SU = ^(abfi\H \aba){a'b'p\H \a'b'u) . (6) 
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The index k labels the various possible combinations of (a,b,a',b'). In the 
direct product space {/z^}, the matrix A^ k ' can be diagonalized and expanded 
into bilinear forms of its eigenvectors ||, 

4£U ^E^WE^^O^l^)) , (7) 

where the A^ s '(k) are real positive eigenvalues, the corresponding right 
(r) and left (7) eigenvectors, r distinguishes degenerate eigenvectors, and 
N(k) > 1 is a suitably chosen normalization factor related to N aa and 

N a (a+1)- 

The generating function Z carries in the exponent terms of the form 

-5EEV 1/, W ff ol flfe ') il/a$ * • ( § ) 

aft 

Here ^ a b^ is a graded vector. The graded matrix L is defined in Ref. [^]. 
After averaging Z, the quartic terms in the ^'s appearing in the exponent 
involve the matrices A^ . The eigenvalue expansion (|7|) allows us to use the 
Hubbard-Stratonovich (H-S) transformation. We introduce the matrices 

xCH(aba'b';p) (9) 

where p specifies whether a right- or left-hand eigenvector C is involved. The 
terms quartic in the ^/'s are bilinear (rather than quadratic) in the matrices 
A. This fact necessitates the introduction of suitable linear combinations of 
the A's. In the H-S tranformation, each such linear combination generates a 
a-field. After integration over the \l/'s, the ensemble average of Z takes the 
form 

(*)\2 



J d[o-]exp -X>g(^(A;)(o- (fc) ) 



trg tr In N(J) 



(10) 



The sum over k extends over all a-fields, the volume element is denoted by 
d[cr]. Here N( J) is both, a graded matrix and a matrix in Hilbert space with 
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basis vectors \abfi) and b = a, (a + 1). The symbol J denotes the source 
terms 0. For lack of space, we give the structure of N( J) only for the cases 
a = a'b = b' where 



(aafi\N(J)\ aav) = ( E - Aa (aa) J^ 

-£Aa^)C^)(0) -w, v (11) 



and 



(a(a+ 1)//|N( J) | a(a+l)i/) = LE - 

_ A(7 (a(a+i)A ^ _ ^Aa^+^C^l) . (12) 



Non-diagonal elements like ((a — l)afi\N(J)\aai/) contain hopping contribu- 
tions. 

Saddle Points. The factors N(k) in Eq. fllQf) obey iV(fe) 3> 1. Hence, 
we use the saddle-point approximation with u = ||. The saddle-point 
equations can be cast into the form of a Pastur equation, 

(abfi\X\a'b'u) = (13) 

pa 



x (abfi\H \abcr) ( — —) (a'b' p\H \a'b'u) . 

\E — W — X ) (ab\a'b');ap 

The matrix X contains all the cx-matrices which occur in N(J). 

Eq. (|14D implies that the solutions aj® are zero or approximately zero 
except for those with k = (a, (a + I), a, (a + 1)) and k = (a, a, a, a). The 
corresponding fields obey the equations 

= (14) 



and 



<4 aa) = tr[(- —^U . (15) 

E -w - AcTsp 

The diagonal solution for a^ a+1 ^ takes the standard form cr^ a ^ a+1 ^ = (E/{2X)) 



±i Ax(E), with Ai(E) = Jl - (E/(2X)) 2 . Here Ai(E) is proportional to the 
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spectral density p sp i(E) in box (a(a + 1)) and has the shape of a semicircle. 
This fact points to a deficiency of the saddle-point approximation: The two 
electrons in slices a and (a + 1) do not interact. The spectrum of each has 
semicircular shape. The spectrum in box (a(a + 1)) is, thus, the convolution 
of two semicircles and not a semicircle. This shows that we must take into 
account loop corrections to the saddle-point approximation. This is done 
below. The invariance of the saddle-point equation Eq. fll4|) under pseu- 
dounitary graded transformations T implies that the saddle-point manifold 
has the form 

O-Wo+l)) = T (a(a+1)) [ E _ . A ^ L j ^ T ( a(a+ i))^i _ ^ 

2A 

Eq. (|15D for contains the two-body interaction w. Diagonalizing w 

and using a simple geometrical construction shows that Eq. (|l^) possesses 
either N aa real solutions (this applies for E < E\ and for E > E 2 ), or, 
for Ei < E < E 2 , (Naa — 2) real and two complex conjugate solutions 
a(E)±iA (E). Thus, 

a (aa) = T {aa) [ ^ _ iA (E) L } (T W )^ , (17) 

where A (E) > is proportional to the spectral density p sp o(E) in box (aa). 
The essential difference between the saddle-point solutions in Eqs. ([16]) and 
fll7D lies in the difference between Aq(E) and Ai(E), i.e, in the different 
spectral densities p sp o(E) and p sp \(E). The interaction w deforms the semi- 
circle Ai(E). Aside from a possible overall shift of the spectrum, w causes 
Aq(E) to be smaller in the centre of the semicircle, and to extend beyond 
the end points (— 2A, +2A) of Ai(E). For a weak interaction, we calculate 
E^ using a second-order perturbation expansion in w which yields 

E h2 = E ±(2 + U 2 )\ . (18) 

Here E = (l/N aa )tr(w) and U 2 = {(l/N aa )tr(w 2 ) - [(l/iV aa )tr( W )] 2 }/A 2 . 
This shows shift and widening of the spectrum. If we identify the matrix 
elements h[f with those of the impurity potential Vi mp , it is easy to see that 
a qualitative change of the spectral density occurs whenever 

(aap\w 2 \aap) av > {ai\{yi mp ) 2 \ai) w (19) 
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where the matrix elements are averaged over \i and i, respectively. This is 
the announced qualitative criterion. 

Localization. To calculate localization properties, we expand the effective 
Lagrangian in the exponent of Eq. (pi]) in powers of v , keeping terms up 
to second order. We put J = and use that for k = (a, a, a, a) and k = 
(a, a + 1, a, a + 1), the ct^'s are solutions of the saddle-point Eqs. ( fH]) and 
(|T5D. We also use Eqs. (|16|) and ([D]). The terms quadratic in the cr^'s 
vanish, and the matrix N yields 

+ ( W /A) 2 (3/ 2 /4)AoA 1 |:trg(r^)L(r^))- 1 
xT 0'+i) L(T (;+i))-i^ (20) 

where the 2K — 1 boxes (aa) and a(a + 1) of the two-electron problem have 
been mapped onto K' = 2K — 1 slices by putting j = 2a — 1 for the boxes (aa) 
and j = 2a for the boxes (a(a + 1). The expression (|20|) has exactly the form 
of the non-linear sigma model derived for the transport of a single electron 
through a disordered wire in Ref. |J. Since Ao and Ai are proportional 
to the spectral densities p sp o(E) and p sp i(E) in boxes (aa) and (a(a + 1)), 
respectively, as determined by the saddle-point condition, Eq. ( |20"D suggests 
defining the saddle-point approximation T sp to the transport coefficient by 

T sp = 2tt p sp0 (E) v 2 Pspl (E) . (21) 

The localization length is proportional to T sp . We conclude that the ratio 
of the localization lengths ((w ^ 0) for non-vanishing two-body interaction 
and ((w = 0) for vanishing two-body interaction is given by 

C(w ± 0) = p sp0 (E,w^0) 

C(w = 0) p sp0 (E,w = 0) • 1 ' 

This result shows that and how the two-body interaction modifies the local- 
ization length. Depending upon the energy considered and on the sizes of the 
shift and of the widening of the spectrum in box (aa) compared to the unper- 
turbed spectrum in box (a(a + 1)), the localization length may decrease or 
increase. In general, we expect an increase near the edge of the unperturbed 
spectrum, and a decrease in the center. 
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Loop Corrections. The result Eq. ( |22| ) is obtained from the saddle-point 
approximation, and it was remarked earlier that this approximation is de- 
ficient. Therefore, we now study the loop corrections to the saddle-point 
solution. We do so by expanding the cr's around the saddle-point solutions 
ai^\ We recall that some of these vanish. We write 

= ffg) + T^Sa^iT^)' 1 ■ (23) 

We expand Z in powers of the Jerk's, keeping in the exponent only the 
quadratic terms which originate from the corresponding terms in Eq. (|10|). 
We need to address the form of the source terms J contained in N(J). The 
expression (|5]) is reproduced if the source term has two non-zero matrix 
elements, (a(a+l)n\ J\(a+n)(a+n+l)u) and {{a+n){a+n+l)v\ J\a(a+l)n). 
The loop expansion generates contributions in which, besides the fork's, 
powers of J and v appear to all orders. We clearly need consider only terms 
of second order in J. We first consider such terms which are of zeroth order 
in v. Inspection shows that, for n ^> 1 and to arbitrary order in the 5(7"'s, 
the result Eq. (|22|) remains valid after Wick contraction of the <W fc )'s. Terms 
of order v 2 come in two classes. There are corrections of order v 2 which either 
appear under the same trace (s) as the source terms to begin with, or which 
become connected to the source terms through Wick contractions of the 5<r's, 
and there are other such corrections for which this is not the case. Formally, 
the latter can be calculated by putting the source terms equal to zero. We 
can show that after Wick contraction, every one of these latter terms takes 
the form 

K ' f 
j=i ^ 

xT (i+i) L ( T (i+i))-i^ ; (24) 

with c a constant which depends on the particular term considered. Re- 
exponentiation of this result yields a renormalization of the transport coef- 
ficient T sp . This is a very satisfying result. Although we are not able to 
calculate the values of the coefficients c analytically, we can show that each 
of these coefficients is determined entirely by two neighboring boxes, say (aa) 
and (a (a + 1)). Any such pair of boxes yields the same result. Thus, our an- 
alytical approach has reduced the calculation of the renormalized transport 
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coefficient to a problem involving only two neighboring boxes, and not the 
full Hamiltonian H . To understand the physical origin and meaning of the 
renormalization terms in this reduced framework, we use the diagrammatic 
expansion of Ref. || which is akin to a diagrammatic impurity perturbation 
expansion. One can see heuristically that the correction terms add up to yield 
the ensemble average of the product of the physical level densities po(E) and 
Pi(E) in the two boxes. Thus, the sum of the saddle-point value and of all 
the loop corrections yields for the renormalized transport coefficient T the 
expression 

7» = 2vr v 2 p {E,w) Pl {E) . (25) 

The ensemble average in Eq. (|25| ) can be determined without difficulty by 
numerical simulation for every given two-body interaction w. We must finally 
take the continuum limit by letting the longitudinal length d of each slice go 
to zero and the number K of slices go to infinity. In this limit we have 

t^^^. (26) 
CO = 0) u T(u; = 0) V ; 

We believe that in a numerical simulation, the ratio of localization lengths 
will be fairly independent of d. 

Summary. We have studied the localization properties of two interacting 
electrons in a quasi one-dimensional wire. Both electrons move in a disorder 
potential and interact when they are close to each other. Hopping allows the 
electrons to move along the wire. 



Our central result is contained in Eqs. (|25|j26|) . The two-body interac- 
tion affects the localization length because it alters the level density in the 
boxes containing two electrons. Depending on the energy chosen and on the 
sign of the two-body interaction, the localization length may decrease or in- 
crease. To the best of our knowledge, this is the first time that an analytical 
treatment of both, disorder and interaction, has been possible. 

The supersymmetry method does not apply to one-dimensional systems. 
Hence it is difficult to compare our results with previous numerical work. The 
exception is the work of Imry who derives an expression (his Eq. (4)) for 
the change of the localization length which applies also in higher dimensions 
and is proportional to (U/B) 2 , with U a typical two-body matrix element 
and B the bandwidth. This expression differs from our expression fl25|) in 
which we retain the structure typical for the Thouless block scaling argument 
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with hopping. We predict a change of localization length which depends upon 
all moments tr(w), tr(w 2 ), ... of the two-body interaction while Imry's result 
involves only his U 2 and is, thus, independent of the sign of w. 



References 

[1] D. L. Shepelyansky, Phys. Rev. Lett. 73 (1994) 2607; see also O. N. 
Dorokhov, Zh. Exp. Teor. Fiz. 98 (1990) 646 [ Sov. Phys. JETP 71 
(1990) 360 ] 

[2] Y. Imry, Europhys. Lett. 30 (1995) 405 

[3] P. H. Song and F. von Oppen, Phys. Rev. 59 (1999) 46 

[4] L. Benet, T. Rupp, and H. A. Weidenmiiller, Phys. Rev. Lett. 87 (2001) 
010601-1; Ann. Phys. (N.Y.) 292 (2001) 67 

[5] S. Iida, H. A. Weidenmiiller, and J. Zuk, Ann. Phys. (N.Y.) 200 (1990) 
219 

[6] K. B. Efetov, Adv. in Physics 32 (1983) 53 

[7] J. J. M. Verbaarschot, M. R. Zirnbauer, and H. A. Weidenmiiller, Phys. 
Rep. 129 (1985) 367 

[8] D. Agassi, G. Mantzouranis, and H. A. Weidenmiiller, Phys. Rep. 22 
(1975) 145 



10 



